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بهینه‌سازی همزمان شکل پوسته و پیکره سخت کننده با روش مجانب‌های پویا* 


۳ (0) ۱) 


بهروز حسنی حسین قاسم نژاد مقری سید مهدی توکلی 


چکیده در اين مقاله بهینه‌سازی شکل و پیکره سخت کننده پوسته‌ها به صورت همزمان انجام شده است. هدف به دس ت آوردن شکل بهینه 
پوسته همراه با چیلمان و نحوه قرارگیری بهینه لایه‌های سخ تکننده در دو سوی سطح سازه پوسته‌ای است. به عبارت دیگر با استفاده از ایسن 
روش همزمان با بهینه‌سازی شکل صفحه میانی پوسته, پیکره سخ تکننده‌های سطحی نیز بهینه شده است. روش حل مسله برمبنای گسسته‌سازی 
اجزا محدود است. مصالح هر جزء توسط مواد ساعنگی با حفره‌های ریزمفیلس مستطیلی شبیه‌سازی شده‌اند. چگالی مصالح هر جز». به صورت 
تابعی از پارامترهای هندسی حفره‌های ریزمفیاس بیان شده و به عنوان متغیر طراحی در فرایند بهینه‌سازی پیکره مورد استفاده قرا رگرفته‌اند. هندسه 
پوسته با استفاده از سطوح نربز تولید شده و تقاط کنترلی تولیدکننده سطح به‌عنوان متغیر طراحی در قسمت بهینه‌سازی شکل استفاده شده‌اند. 
مسئله بهینه‌سازی در هر بخش با استفاده از روش مجانب‌های پویا حل شده است. به منظور نشان دادن کارایی روش, چند مال ارائه شده و نتایج 


مورد بررسی قرا رگرفته‌اند. 


واژه‌های کلیدی بهینه‌سازی پیکره. بهینه‌سازی شکل پوسته سخت کننده. 
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۳۸ 


مقد مه 

بهینه‌سازی سازه‌ای در دهه‌های اخیر پیشرفت‌های زیادی 
کرده و کاربردهای موفقیت‌آمیزی در مسائل مهندسی 
داشته است. در بهینه‌سازی شکل. با مفروض بودن 
بارگذاری و شرایط مرزی» هدف دست‌یابی به بهترین 
هندسه مرزهای دامنه طراحی است. از جمله مراجع در 
این زمینه می توان به زینکویچ (2160160160102) و کمپیل 
(ااعطمصعت) [1] هافکا (11211) و جاندلی (زطهدم) 
[2] و دینگ (ه۳) [3] و همچنین اولهف (01۳08) و 
رزوانی (07۷201) [4] اشاره نمود. 

در بهینه‌سازی شکل به صورت متعارف آن فرض 
بر آن است که پیکره (۱0001081) دامنه مسئله ثابت است. 
اما در حالت کلی روش‌های بهینه‌سازی که فادر به ایجاد 
تغییر در تعداده موقعیت و نحوه ارتباط عناصر سازه‌ای یا 
به طور خلاصه پیکره سازه باشند. مورد نیاز هستند. در 
سال‌های اخیر روش‌های متنوعی برای این نوع 
بهینه‌سازی معرفی شده‌اند و در حال حاضر روش‌های 
زیادی برای حل مسئله پیکره وجود دارد. در تحقیقات 
اولیه که توسط بندسو (8000500) و کیکسوچی 
(ندهدان) [5] انجام روش حل بر مبنای استفاده از 
ریزسازه‌ها (معسهنااه ۷:6۲0) و به دست آوردن 
خواص مکانیکی مصالح با استفاده از تشوری همگن 
سازی (۲10000860122100) بوده است. یک روش 
جایگزین برای آن» استفاده از مصالح ساختگی همسانگرد 
بااعمال تاوان با ۹۲/۴ ( ۷۲26۵۲121 15000010 50110 
حمناهعنآهمع۳ ط1ز) است [6]. در این روش ماده مورد 
استفاده همسان فرض می‌شود و چگالی به عنوان تابعی از 
مختصات مصالح مقداری بین صفر و یک اختیار می‌کند 
طوری که چگالی یک به معنای وجود و چگالی صفر به 
معنای عدم وجود مصالح بیان می‌شود. علاوه بر ایین 
روش‌ها. روش‌های دیگری نیز برای بهینه‌سازی پیکره 
وجود دارد. از آن میان می توان به روش‌های تکاملی 
(ممنلمء‌نستامه . بعطه‌نابآ۳۵) (ژی و استیون [7]) 


روش کلونی مورچه‌ها (۰01057 ۵۸۲) [8]» روش سازه 
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مبنا (متستعتاه فصم0) [9]؛ و دیگر روش‌ها اشاره 
نمود. 

به طور معمول در بهینه‌سازی پیکره نحوه توزیع 
مصالح در یک دامنه ثابت جستجو می‌گردد و سازه به 
دست آمده به عنوان سازه پایه جهت بهینه‌سازی شکل 
مورد استفاده قرار می‌گیرد (اولهوف و همکاران [10]؛ لین 
ای ار سا سکره سکف ککارم اوه 
پوسته‌ای در حالتی که شکل صفحه میانی آن (سازه پایه) 
ثابت است لزوماً به معنی رسیدن به بهترین پاسخ یا حتی 
یک پاسخ مناسب نیست. زیرا در نظر گرفتن شکل 
جدیدی برای سازه پایه ممکن است منجر به پاسخ بهتری 
برای توزیع مصالح شود [12]. آنسولا (۵۵5014) و 
همکاران [12] روشی را معرفی کردند که در آن از یک 
فرایند ترکیبی برای بهینه‌سازی همزمان شکل و پیکره 
استفاده می‌شود. در این روش مسئله ترکیب شده به دو 
زیرمسئله شکل و پیکره تقسیم می‌شود. از این‌رو از یک 
روند تکراری دو مرحله‌ای برای بهینه‌سازی استفاده 
می‌شود؛ به طوری که در مرحله اول با انجام یک گام 
بهینه‌سازی شکل. هندسه پوسته بهبود می‌یابد و در مرحله 
دوم توزیع بهینه لایه سخت کننده روی پوسته با استفاده 
از روش معیار بهینگی انجام هی گراخ این کار تا همگرا 
شدن به یک پاسخ بهینه ادامه پیدا می‌کند. 

ایده اصلی فن توزیع مصالح برای مسائل 
بهینه‌سازی پیکره تبدیل مسئله به مسئله بهینه‌سازی اندازه 
برای چگالی در محدوده ثاببت است. مشکل عمده در 
ینجا پیدا کردن الگوریتم‌های بهینه‌سازی است که قادر به 
حل مسئله با تعداد متغیر طراحی زیاد باشد. بدین منظور 
ز دهه شصت میلادی روش‌های معیار بهینگی در حسل 
مسائل بهینه‌سازی سازه‌ای مورد استفاده قرار گرفشه‌اند. 


ین روش‌ها. روش‌های غیر مستفیم بهینه‌سازی هستند و 


هدف را بهینه می‌کنند» تلاش می‌کنند تادسته‌ای از 
معیارهای مربوط به رفتار سازه را اقناع نمایند. در 


سال‌های اخیر استفاده از روش‌های برنامه‌ریزی ریاضی 
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بهروز حسنی- حسین قاسم نژاد مفری- سید مهدی توکلی 


در مسائل بهینه‌سازی پیکره همانند مسائل بهینه‌سازی 
لاخ و شک یشان وله تم رفن آش هویش هت 
مسئله طراحی به عنوان مسئله بهینه‌سازی بر حسب 
متغیرهای طراحی در نظر گرفته شده و توابع هدف و قید 
بر تیان طقف قاط ای مان ی تفن فاد 
از مقادیر مشتقات توابع هدف و قید بر حسب متغیرهای 
طراجیها ییاز ی حل فی‌سود: 

در اين مقاله مسئله بهینه‌سازی همزمان شکل و 
پیکره لایه سخت کننده سازه پوسته‌ای به صورت مسئله 
حداکثر سازی سختی سازه با استفاده از حداقل‌سازی 
کارخارجی انجام شده روی سازه در نظر گرفته شده 
است. شیوه کلی حل مسئله بر اساس روش دو مرحله‌ای 
پیشنهادی آنسولا و همکاران است؛ طوری که در هر 
تکرار بهینه‌سازی در مرحله اول شکل سطح پوسته و در 
مرحله دوم توزیع لایه سخت کننده آن بهینه می‌شوند. 
این مراحل تا همگرایی پاسخ متناوبا تکرار می‌شوند. برای 
بهینه‌سازی از روش ریاضی مجانب‌های پویال[ ۱۳] استفاده 
شده است. این الگوریتم نشان داده است که تطبیق‌پذیر 
بوده و برای مسائل با تعداد متغیرهای طراحی زیاد 
مایت 2 


مدل مصالح در بهینه‌سازی پیکره 

معرفی یک تابع برای مشخص کردن چگالی مواد در 
جسم با در نظر گرفتن مواد مرکب مشتمل بر تعداد 
نامحدودی از سوراخ‌های ریز که به صورت متناوب در 
این جسم پخش شده‌اند. مسئله بهینه‌سازی پیکره سازه‌ها 
را به یک مسئله بهینه‌سازی آندازه تبدیل می‌کند و 
پیچیدگی این مسائل تا حدودی برطرف می‌شود. در واقع 
اتفاهه ان یهگا گیری ی رانا ریرسازهاق 
متناوب (وعستاعب۹ معنا عت۳۵۲۵)؛ طبیعت بودن با 
نبودن مواد در پیکره جسم را از مقیاس ماکروسکپی به 
مقیاس میکروسکپی رهنمون نموده و باعث پایداری حل 
می گردد. 
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یکی از راه‌هایی که برای مدل کردن مواد در مسائل 
گسسته بهینه‌سازی پیکره وجود دارد استفاده از مدل 
مصالح ساختگی همسانگرد با اعمال تاوان است. در این 
مدل بودن يا نبودن مواد با استفاده از تابع تقریب مناسبی 
که تابع چگالی مصالح ساعتگی نامیده می‌شود. مشخص 
می‌گردد. این تابع مقادیر بین صفر و یک دارد به طوری 
که تابع چگالی صفر به معنی نبود مصالح و مقدار یک به 
معنی وجود مصالح می‌باشد. با استفاده از این مدل 
می‌توان خواص مکانیکی ماکروسکپی مصالح را بدست 
آورد. 

در انتخاب ریزسازه‌ها. شکل آنها یکی از مهمترین 
مسائلی است که باید به آن توجه کرد. شکل انتخاب شده 
باید به گونه‌ای باشد که چگالی مواد در ریزسازه بتواند 
کل مقدار ۰ تا ۱ را پوشش دهد. از طرف دیگر شکل 
حفره باید با کمترین تعداد پارامتر تعریف شود تامیزان 
متغیرهای طراحی در مسئله بهینه‌سازی به حداقل برسد. 
سلول‌های مربعی پا حفره‌های مستطیل شکل در مرکز آنها 
ساده‌ترین شکل برای این منظور می‌باشند (شکل (۱)). 


شکل ۱ سلول‌های مربعی با حفره‌های مستطیلی 


اگر بعد سلول مربعی با 8 که مقدار مثبت و بسیار 
کوچک است نشان داده شود آنگاه اندازه سوراخ با 82 و 
0 مشخص می‌گردد. شکل (۲) سلول واحد را در 
دستگاه مختصات میکروسکپی نشان می‌دهد. با استفاده از 
این مدل. سطح اشغال شده توسط مواد جامد به‌صورت 


زیر محاسبه می‌شود. 
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شکل ۲ سلول واحد با حفره مستطیلی در مختصات میکروسکپی 
)۱( 0 (01-0 ی < ,62 


که در ایین رابطه 41 0>8 و 04961 و ۵ فضای 
طراضی و 6 قسصت بحابل فشای ای ۷ فان 
می‌دهد. با توجه به نکات فوق هر نقطه 2 > » دارای 
فقتاویر الق تروش تاهت کته معفی هی بط رای وس 
بهینه‌سازی پیکره هستند. 
)۲ (0 ع<ه , (20< 2 
در عمل این توابع به‌وسیله توابع ثابتی در هر جزء 
از فضای طراحی گسسته سازی‌شده تقریب زده می‌شوند 
و بنابراین ابعاد یک ریزسازه در هر جزء ثابت فرض 
می‌شود. در نتیجه ماتریس کشسانی همگن شده نیز برای 
هر جزء ثابت خواهد بود. بنابراین در فضای دوبعدی اگر 
دامنه به لا جزء محدود تقسیم شود 21 متغیر طراحصی 
در مسئله بهینه‌سازی پیکره سازه وجود دارد. تابع چگالی 
ساختگی ی برای ساختن ریزسازه‌های ساختگی با 
پارامترهای هندسی در ارتباط است. برای ساختن جسم 
سلولی شامل سلولهای واحد با حفره‌های مستطیل شکل 


۰ (06) به‌صورت زیر در نظر گرفته می‌شود: 
(۳( 000 21-200 (500 


در این رابطه 1 > (*)8 > 0 و 62 > * می‌باشند. به 
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ر ۵( - ,62 


همچنین چگالی و ماتریس کشسانی به صورت زیر 


خواهند بود: 


20۲ - ()م 
(20 - 000 


2 
"م و ") به ترتیب چگالی و ماتریس کشسانی 
قسنت جامد همگن می‌بافنند. درز این رابطه و غتریب 
خر یه و یفن کین ۱ رمع شم ۲و۹ کاس اشتا: 
همان‌طور که اشاره شد در مدل به‌کار رفته فرض می‌شود 
که مصالح حاوی تعداد زیادی سلول‌های مربع شکل با 
حفره‌های مستطیلی بوده و پارامترهای (200 و (9)۶ در 
هر جزء محدود ثابت هستند. حال مدل سلولی معرفی 
می‌شود که از آن در فرایند بهینه‌سازی پیکره لاه سخت 
کته مرا ساه‌های :یر شعهای اشفا ده اه قنل یا 
استفاده از نظریه مونتاژ لایه‌ها (۳۱65 عصناطصعه‌عه) [15] 
مدل سلولی شکل (۳) را می‌توان مورد استفاده قرار داد. 
در این مدل لایه همسانگرد میانی معرف سازه اصلی و 
لایه‌های بالایی و پایینی تشکیل دهنده لایه سخت کننده 
هستند. برای لایه سخت کننده مدل میکروسلولی با 
مزمای شمان سورد بستانه فر کرو مان گر 
که قبلا اشاره شد ابعاد حفره‌ها و جهت قرارگیری آنها 
پارامترهای سلولی خواهند بود. باید توجه داشت که در 
اتف فرصم سرانت استه کته یبا وس ی انساد 
حفره‌ه‌ای میکروس‌کوپی (0,۵۸) در حین فرایند 
بهینه‌سازی» و) ضخامت سازه اصلی و ,4 ضخامت لایه 
سخت کننده ثابت باقی می‌مانند. در نتیجه اندازه» محل و 
نحوه قرارگیری لایه سخت کننده (با یک ارتفاع ثابت) به 
وسیله یک الگوریتم بهینه‌سازی پیکره. مشخص خواهند 
شد. به خاطر حفظ تقارن لابه سخت کننده در هر دو 


طرف سازه در نظر گرفته می‌شود. 
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شکل ۳ سلول پایه به همراه لایه سخت کننده 


حجم اشغال شده توسط لابه‌های سخت کننده را 


می‌توان از رابطه زیر به‌دست آورد: 
)‌ 0(۵ه-1) مارا ۷ 


که 9۳ محدوده لابه میانی از سازه اصلی است. 


تعریف هندسه مسئله 
تعریف هندسه است. در تحقیق حاضر از سطوح نربز 
(مصتام6-ظ امعمتم؟ صصمئتصنا ممل :۱/886 برای 
تولید سطح پوسته استفاده می‌شود. این کار امکان ایجاد 
شکل‌های پیوسته با استفاده از تعداد محدودی متغیر 
طراحی را فراهم می‌سازد. همچنین این امکان را می‌دهد 
تا بعد از بهبود متغیرهای طراحی در روند بهینه‌سازی. 
شکل جدید پوسته با حفظ پیوستگی و نیز مدل اجزای 
محدود جدید تولید شوند. در این بخش نحوه تعریف 
سطوح هندسی با استفاده از فن نربز مختصرا شرح داده 


ی زگ 
نوابع پایه بی اسپلاین 


اوه ۱ مجموعه‌ای پیاپی 51 غیر نزولی از 
اعداد حقیقی است» یعنی در آن به ازای 1- 0,....90 - 1 


سال بیست و سوم شماره دی ۱۳٩۱‏ 


۱ 


داریم لا 2 ,یلا . به بل گره و به مجموعه لآ بردار گره‌ای 
از درجه ۲ (مرتبه 1+) که با (۲) ,۷( نشان داده 


می‌شود به صورت زير تعریف می‌شود [16]: 


0 برزلا > اک ولا 1 1 ۱ ۳ 
00 0 1,0 
1 :11-0 
اب 1( ۸ 
و 9 سا وا تاصا 


برخی از خواضن. توانع پایه‌ای: بسن اسپلاین به 

صورت زیر می‌باشند: 

۱- برای همه مقادیر 0 بو 1> ۷ > 0 داریم: 20 (۷)م :ل. 

7- جمع مقادیر این توابع در 05۱51 برابر واحد 
است: پمتی )و۷ 

۳ در صورتی که کار بازه تاثیر (,بمبزلا,:"] 
باشد. 0 (0) م2 . 

6- بجز در حالت 0<. (0), :لا تنها دارای یک نقطه 
ماکزیمم می‌باشد. 

قهر نک کرهاز تفای گرضی که یار #زان تام ناشن 
(), ,۰۷ 0-1 بار بصورت پیوسته مشتق‌پذیر 
پاش شیاین ا‌انقن ده یات اف رات 
پیوستگی و افزایش تعداد تکرار هر گره سبب کاهش 
پپوستگی و اهل شل, 

1- مشتقات توابع پایه بی اسپلاین را می‌توان به صورت 


زیر نوشت: 


)۸ () بلط ۳ (۲) ولا 1 


زبنلا - تبمبزلا نا وبز 


شکل ( توابع پایه‌ای بی- اسپلاین از درجات ۰ 
۱ و ۲ را نشان می‌دهد. 
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1 

505 "2 

۳ 

۳ خسبام وتو و و او ون و از ال ویو و 
05 ۳ 025 8 1 015 5 
۷ ‌" 
1 

3 
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ِ 
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۳4 
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7 ۷-0 ۳ 
نج 505 ۹ 


‌ ۱ ۱ ۱ 
یی سس اس 
9 7 12 1 


یر 
025 0 


(۱۷۱۵)۷ 
هك 
__ ه 
ه‌ِ ۱ 
ط 
‌ِ : 
۷ 


شکل ۶ توابع پایه بی اسپلاین از چپ به راست از درجه صفر تا دو 


سطوح نربز 

از توابع پایه بی اسپلاین می‌توان برای تولید منحنی‌هاو 
سطوح هندسی پیچبده استفاده نمود. به سطوحی که با 
استفاده از توابع پایه بی اسپلاین به صورت کسری و 
اعمال ضرایب وزن تولید می‌شوند. سطوح نربز گفته 
می‌شود [16]. 

یک سطح نربز که در جهت "از درجه و در 
جهت ۷ از درجه و باشد. را می‌توان به وسیله توابع 
کسری قطعه‌ای به صورت زير تعریف کرد: 


گر ۷(۷) رل (0) وی ۱ 2 
20 


5)0,:۷( :0 > ۲,۷ >1 . )0( 


زر (۷) و لا (۲) ,رل( 


0-ز 1-0 


در عبارت فوق [,,۳) شبکه نقاط کنترلی است که 
در دو جهت تعریف شده است؛ همچنین ارت وزن‌ها. 
((«) ,,) و ((") ,رل توببع پایه‌ای بسی- اسپلاین 
هستند که بر روی بردارهای گره‌ای رابطه (۱۰) تعریف 
شده‌اند. در این روابط تساویهای 0+1 +90 < و و 
0+1 +7 < ۲ برقرار می‌باشد. با تعریف توابع پایه‌ای 
نسبی قطعه‌ای به صورت رابطه (۱۱). می‌توان رابطه )٩(‏ 
تاه موی (۱۷) رک تیه گر قسب یا اسهاوی 3۱ 


رابطه (۱۲) می‌توان هرگونه سطح و رویه پیچیده‌ای را 


نشریه مهندسی عمران دانشگاه فردوسی مشهد 


مدل نمود. 
او وف لا یلا ولا اکتا 
2+1 ۱+ 
۱ ۰ ویو :۷ دربو ۰0۰۷ 0 
آب+و + 
۰ 
۷ از 
6۱ انریا ( )ود << (۷:) ر .18 
62 () شور 
1-0 1-0 
٩)0,۷( < ۱ 3 ):۷( ۲ ۱۲(‏ 
0-ز 0-ز 


طبیعت پارامتریک مدل طراحی امکان وارسی شکل 
سازه را از طریق بهبود محل قرارگیری نقاط کنترلی ایجاد 
می‌کند. به منظور کاهش تعداد متغیرهای طراحصی. 
می‌توان تنها تعداد محدودی از نقاط کنترلی را به عنوان 
متغیرهای طراحی در نظر گرفت. با حرکت نقاط کنترلی 
در جهت‌های مشخص. شکل‌های متفاوتی به دست 
می‌آید. محل‌های قرارگیری نقاط کنترلی به عنوان 
متغیرطراحی در نظر گرفته می‌شوند. متغیرهای طراحی 
بهینه‌سازی شکل, مستقل از هم هستند و هر مجموعه از 
آنها منجر به ایجاد طراحی متفاوتی از شکل خواهد 


شد. 
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رابطه سازی مسئله بهینه‌سازی 
یک مسئله کلی کشسانی خطی تحت اثر نیروهای حجمی 
به‌ کار برده شده در دامنه 62 و نیروهای سطحی در ,1 
مفروض است. سطح دامنه 62 دارای مرز 1 شامل ,1 که 
تغییرمکانها در آن تعریف شده‌اند و ,۲ که بارهای سطحی 
در آن به کار برده شده‌اند می‌باشد. همچنین فرض 


می‌شود که: 
دواد ] 
20 1 [ 


۱۳ 
با به‌کارگیری روش تغییرمکان مجازی, معادلات 
تعادل با مساوی قرار دادن کار مجازی داخلی و خارجی 
به دست آورده می‌شوند. با فرض اینکه ۷ میدان 
تغییرمکان که تعادل سازه الاستیک را تعریف می‌کند باشد 
و ۷ میدان تغییرمکان مجازی مجاز سینماتیکی باشد 
,]9 2-0 ۷ 20704 (6) )2 ۷ ۷] > ۷ ۲0عط ۰۷ ۷ 

(۱ 
پس برای یک سازه کشسان با مرز ثاببت و1 


 6۵(‏ 17۷0۲ ]+۱۵۵ ‌- (۱(6۵۵) " أ 


با بهره‌گیری از تحلیل حساب تغییراتی و تابع نماها 
و استفاده از فرم انرژی دوخطی برای کار داخلی و فرم 
بارخطی برای کار خارجی رابطه (۱۵) را می‌توان به 


2)0,۷( 2 /)۷( :, / ۷ (۱۹) 

که در این رابطه: 
۳( 6 () قأ < (۵)0,۷ 
0۸۳ ۱۷ ]| +1۵ - 10 
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1۳ 


هدف از مسئله بهینه‌سازی سازه‌ها که در این مقاله 
رف مت فان کر امس تاه 
ممکن با استفاده از یک مقدار مصالح معین است. مینیمم 
بودن مقدار کار خارجی با میدان تغییرمکان حقیقی و یا 
مینممم بودن (۰/6 ماکزیمم بودن سختی عمومی یک 
سازه را فراهم می‌کند. بنابراین مسئله بهینه‌سازی سازه‌ها با 


قرار دادن (0۲)/ به عنوان تابع هدف به صورت زیر ساخته 


ی نود 
 200(‏ ۱۷۲1919176 
۷۷ (۷// 2 (200,:۷ .5 
عصمتام‌تتاوع؟ طعنوع0 ...۰ 2۳00 
)۱۹ 


که در این رابطه محدودیت‌های طراحی» (۸)0,۷ و 
و متغیرهای طراحی با توجه به نوع مسئله بهینه‌سازی 
(شکل یا پیکره) تعیین می‌شوند. 

دو مرحله‌ای پیشنهادی آنسولا و همکاران [12] استفاده 
می‌شود. در مرحله اول. بهینه‌سازی شکل انجام می‌شود و 
متغیرهای طراحی بهینه شده که شکل سطح میانی پوسته 
را تعریف می‌کنند به‌دست می‌آیند. این متغیرهای طراحی 
همان نقاط کنترلی تولید کننده سطوح نربز هستند. 
محدودیت طراحی که در این مرحله وجود دار حدود 
تغییرات نقاط کنترلی است که توسط طراح تعیین 
می‌شود. مسئله بهینه‌سازی در این حالت شامل 9 متغیر 


۱۷۲9126  0( 
5. 200:۷( 2 0۷/۳ ( ۷ 
220 رراحک رک > ف‎ 
۲۰( 


کین ان 2 محل نقطه کنترلی و 9 و 0 حدود 
در مرحله دوم چگالی‌های مصالح که برحسب 
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در هر جزء ثابت‌اند. برای توزیع ماده تشکیل دهنده 
سخت کننده پوسته بهینه می‌شوند. قید طراحی که در این 
قسفت. اعمال می‌شود. مقدار مصالح مورد استفاده برای 
شود. بدین ترتیب مسئله بهینه‌سازی شامل 1 متغیر 
طراحی (9 تعداد جزء‌های محدود است) و یک قید به 
صورت زير خواهد بود: 


۱۷۲1019176  00( 
8 


5.1 200:۷( < //۷( ۷ 


۳ 1 
۰ > 062 زم | .ده 21010 
(۵ 
1-1 
هرسیاحف > رو > مرن > 0 


1 
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که 6 تک بل بحج: مط الم اما نی با شتا 
همچنین ,زمر حد پایین چگالی است که برای 
جلوگیری از تکینی ماتریس سختی اعمال می‌شود 
که در عمل مقدار آن 20.001 زیم قرار داده می‌شود 
[14]. 

روند حل پیاپی مسائل شکل و پیکره تا همگرایی 
به یک پاسخ درست ادامه می‌یابد. به منظور جلوگیری از 
تغیبرات زیاد در طراحی بهتر است تعداد تکرارها در 
حلقه‌های داخلی کم باشد [12]. برای مثال‌های حل شده 
در این مقالهء تعداد تکرارها در هر حلقه برابر یک در نظر 
گرفته شده است. شکل (۵) نمودار جریان بهینه‌سازی 


همزمان را نشان می‌دهد. 


بت سس سس سس سس جسست بت ستت. بت تست سس سس بت مس م۲ 


شکل ۵ نمودار جریان بهینه‌سازی همزمان شکل و پیکره 
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روش مجانب‌های پویا 
روش مجانب‌های پویا و روش مادر آن حطی سازی کوژ 


]17[ (۷۰: 0۵0 ].]۵۵۲1721100( 


الگوریتم‌های برنامه‌ریزی ریاضی هستند که برای مسائل 
طراحی پیکره مناسب‌اند. این روش‌ها همانند روش‌های 
حطی سازی یر 


8 و سهموی سازی پیاپی ( :90۳ 


پیاپی ( تهعصنا امتلعدوعک ٩81۳:‏ 


ومتصمصموم ناملعم امتامعناوع8) برای حل مسائل 
بهینه‌سازی هموار غیرخطی از دنبال‌ای از زیرمسئله‌های 
تقریبی مشابه استفاده می‌کند. برای مجانب‌های پویاو 
خحطی سازی کوز این زیرمسئله‌ه ا تفکیک‌پذیر 
(960272010) و کوز (600۷96) می‌باشند و بر اساس 
مقادیر مشتقات در نقطه تکرار جاری و برخی نقاط تکرار 
قبلی ساخته می‌شوند. در هر نقطه تکرار این زیرمسئله با 
اوه درون درگتان [19] سای هی رون بانیم 
زیرمسئله به عنوان طرح بعدی در فرایند تکراری مورد 
استفاده قرار می‌گیرد. 


در روش مجانب‌های پویا تقریب تابع ۴ باه 


متغیر حقیقی (,ا,...,) < 2 در اطراف نقطه تکرار 

پا سکن ویو انسیت: 

(۲۲ ی +( ۳0 د ۲00 
( مت ج) ی *( 0 


که مقادیر 1 و :5 بصورت زیر می‌باشند: 


20 9 : (00 رد لا ود 0 < من 2 ۶ 


1 1 


7 
ِ (نا- )-- ره ,0 وج > (0من 2 م 


مقادیر مثبت لا و نا وارسی کننده‌ی محدوده‌ای 
هستند که در آن تقریب 8 قادر به تولید پاسخ‌های 
معقول برای مسئله بهینه‌سازی است (پارامترهای زلاو ,1 
مجانب‌های قائم برای تقریب هستند). این مقادیر در 
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الگوریتم بهینه‌سازی, برای هر یک از توابع مسئله در هر 
تکرار بسته به تاریخچه تکرار به روز می‌شوند. 

بدیهی است که افزايش چگالی. منجر به افزایش 
سختی و در نتیجه کاهش میزان کار خارجی می‌شود. 
بنابراین همواره مشتق کار خارجی برای هریک از 
چگالی‌های هر یک از اجزا منفی خواهد بود. در نتیجه 
تقریب کار خارجی در تکرار گام زیرمسئله به صورت 


(ام) ۳ 2[ کر (*م)/) هنده 


.00 سل و 1 


12 
تا 0۵ ی 6 > ۰( 5 
دز 


(۲۶ 
حل ا ین مسئله با استفاده از روش دوگان امکان 
پذیر است. ابتدا باید تابعی (۲۵6/021) لاگرانژ بآ به 


صورت زیر تشکیل می‌شود: 


زک مر (*م/ - ,1 


۸ )60,-60,( 


اح[ 

۹۳2 
سپس این تابعی نسبت به چگالی‌های :م که 
8ر...یع زر 1گ ز۵ ک ونیم حداقل می‌شود. به خاطر 
کوژ بودن و تفکیک پذیر بودن سل مسئله به سادگی به 
صورت جزء به جزء حبل می‌شود. گام دوم در روش 
دوگان حداکثرسازی تابعی به دست آمده نسبت به ۸ 
است که این مرحله منجر می‌شود به تنظیم مقدار ۸ 
طوری که الگوی بهبود منجر به بدست آوردن مقادیر 
شتبری کته قیل یی وا تاسن کته می کوان 
مقادیر زا را به صورت هوشمندانه‌تن طوری انتخاب 


نمود که سرعت همگرایی را افزایش دهد. 
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حل مثال‌های نمونه 
در این بخش به حل مثال‌هایی از روش ترکیبی فوق و 
مشاهده نتایج آن پرداخته می‌شود. در شکل‌های این 
بخش رنگ سفید بیانگر سازه اولیه و نواحی سیاه شده 
نقان دفنده تراحی تقویت شدهس پاش ساز‌ان: :یه 
عنوان اولین مثال مورد بررسی قرار می‌گیرد پوسته‌ای به 
شکل نیم‌استوانه است که چهار تکیه‌گاه ساده در 
گوشه‌های آن قرار دارند. از یک زیر برنامه تولید سطح با 
استفاده از منحنی‌های نربز برای ایجاد هندسه سازه 
استفاده می‌شود. بدین منظور از دوازده نقطه کنترلی برای 
تولید شکل استفاده شده است که مختصات قائم محل 
قرارگیری آنها متغیر طراحی در مسئله بهینه‌سازی شکل 
خواهند بود. شکل (۱) نحوه قرارگیری نقاط کنترلی و 
هندسه اولیه سازه. به همراه شبکه‌بندی آن را نشان 


می‌دهد. 


شکل 1 سازه استوانه‌ای به همراه نقاط کنترلی تولید کننده 


سطح آن 


برای تحلیل به روش اجزا محدود سازه به ۸۸۶ 
جزء محدود تفسیم شده است. ضخامت سازه پوسته‌ای و 
لایه سخت کننده آن هریک 510757 و خصوصیات 
مک‌انیکی مصالح که / 18 21*10۴ < ۴ و 9-03 
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پهینه سازی همزمان شکل پوسته و پیکره سخت کننده 2 


مختصات قائم نقاط کنترلی تولیدکننده سطح هستند. 
بنابراین در این مرحله دوازده متغیر طراحی به همراه یک 
قید محدوده حرکت برای هریک وجود دارد. در قسمت 
بهینه‌سازی پیکره به تعداد هر جزء یک چگالی وجود 
دارد. در نتیجه ۸۸۶ متغیر طراحی به همراه یک قید 
محدوده تغییرات (بين ,زرم و ۱) برای هرکدام وجود 
دارد. علاوه بر آن یک قید حجمی هم در این مرحله 
لحاظ می‌شود که کنترل‌کننده میزان مصالح مورد استفاده 
جهت توزیع در لایه سخت کننده است که به صورت 
درصدی از مصالح کل دامنه بیان می‌شود. برای مثال‌های 
این مقاله. اين میزان برابر ۶۰ درصد در نظر گرفته شده 
است. همچنین ضریب تاوان برای مدل مصالح ساختگی 
برابر 23 0 درنظر گرفته می‌شود. 

سازه مذکور تحت بارگذاری‌های مختلف قرار 
گرفته است که شامل بار متمرکز 51 وارد بر مرکز بار 
گسترده خطی ثابت در وسط به اندازه 150/0 و بار 
گسترده سطحی *2(/000روی کل سطح آن می‌باشد. 
نتایج بهینه‌سازی همزمان حالات بارگذاری مذکور به 
ترتیب در شکل‌های (۰۷ (۱۰) و (۱۲) قابل مشاهده 
است. همان‌طور که ملاحظه می‌شود برای هر یک از 
حالات بارگذاری علاوه بر تفاوت پیکره لابه سخت 
کننده. شکل متفاوتی نسبت به حالت اولیه به دست 
می‌آید که معرف بهترین شکل و بهترین توزیع ماده در 
لایه سخت کننده می‌باشد. همچنین روند همگرایی پاسخ 
در تکرارها برای حالات بارگذاری مذکور به ترتیب در 
شکز ماع (۸( ۱۱و (۱۳) فان داد فله اش: 

در مثال بعدی سازه کره‌ای شکل مورد بررسی قرار 
می‌گیرد (شکل ۱۶). این سازه در لبه انتهایی خود تحت 
تکیه گاه ساده قرار گرفته است. برای تولید و وارسی 
شکل آن از ٩‏ نقطه کنترلی استفاده شده است. تعداد اجزا 
محدود برای این سازه ۲۰۶۸ در نظر گرفته شده است. 
ضخامت و حصوصیات مکانیکی همانند مثال قبلی لحاظ 
شاه است, ار این فذاد متیر های طرای مسافل شکل 


و پیکره در این مثال به ترتیب برابر ٩‏ و ۲۰۶۸ می‌باشد. 


سال بیست و سوم. شماره دی ۱۳۹۱ 


بهروز حسنی- حسین قاسم نژاد مفری- سید مهدی توکلی 


همانند مثال قبلی حجم مصالح تشکیل دهنده لایه سخت 
کنتله. »8 عرص و ضریب کاواخ براب ۳ در نظر گرفته شده 
است. 

دو حالت بارگذاری برای ایین سازه شامل بار 
متمرکز برابر 516 وارد بر مرکز و بار گسترده 
سطحی *50(۷/0 در نظر گرفته شده است. نتایج حاصل 
به ترتیب در شکل‌های (۱۵) و (۱۷) نمایش داده شده 
است. همچنین در شکل‌های (۱۱) و (۱۸) نمودارهای 


۷ 


روند و آیی خمه ثر تسب نشان داده شده‌اند. جنانجه 


مشاهده می‌شود. برای شکل‌های این مثال» شکل اولیه 
تفه کی فان که ات آیی مسگله یه عناظر 
تقارن شکل هندسی کره و بارگذاری‌های این مثال 
می‌باشد. به عبارت دیگر شکل اولیه کره برای 
با گذازی‌قاق انم لباقت آشا گر سفیتم 


کننده با توجه به نوع بارگذاری‌ها تغییر کرده است. 


ی 
۱ 
۳ 1 
و 
90110 
۹۹ 


٩ تیر‎ 


)( 


شکل ۷ بهینه‌سازی همزمان سازه استوانه‌ای تحت بار متمرکز, (الف) نمای کلی» (ب) و (ج) نمای جانبی, (د) نما از بالا 


7500 
7000 7 
6500 
6000 
5500 5 
6000 7 
4500 5 
4000 


آنرری کرنسی 


)۱۱-6( 


شکل ۸ نمودار همگرایی بهینه‌سازی همزمان سازه استوانه‌ای تحت بار متمرکز 


سال بیست و سوم شماره دی ۱۳۹۱ 
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شکل ۱۰ بهینه‌سازی همزمان سازه استوانه‌ای تحت بار گسترده خطی. (الف) نمای کلی. (ب) نما از بالاء و (ج) نمای جانبی 


7500 
7000 
68500 
6000 
5500 
2000 
1500 
4000 
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انرری کرنسی (۱۱-610) 


شکل ۱ نمودار همگرایی بهینه‌سازی همزمان سازه استوانه‌ای تحت بار گسترده حطی 
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سال بیست و سوم. شماره دی ۱۳٩۱‏ 


شکل ۱۶ سازه 


کره 


ی 


شکز 


فا 
۷ 


۳۹ 


20 


15 


شکل ۱۲ بهینه‌سازی سازه استوانه‌ای تحت بار گستر 


ده 


به روز حسنی- حسین قاسم نژاد مفری- سید مهدی ت وکلی 
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شکل ۱۳ نمودار همگرایی بهینه‌سازی همزمان سازه استوانه‌ای تحت بار گسترده یکنواحت 


۵ 
ی زا 95 9 95 9 و9 ها 
هب با هب با فا با ها ف 
۵ (۶ا ۵ (ا ی شا تا ها 
۰ ۵ ما هه ف8 ۲ ۲ 
(۷۱۵-۱۱) کج کی ۳۲۶7 


یکنواخت. (الف) نما از بالاه (ب) نمای جانبی. و (ج) نمای کلی 


1۹ 
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شکل ۱۵ بهینه سازی سازه کره‌ای تحتباو متمرکو (الف) نما جانبی (ب) نماق از بالا 


11000 
10000 
9000 17 
8000 7 
7000 
6000 
5000 


000 1 1 
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نرزی کرنسی (۱۸-2۳0) 


شکل ۱۸ تمودار همگرایی بهیته‌سازی عمزمان سازه کره‌ای تحت بار عشمرگر 


424 
۳7+ 
۸93 
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شکل ۱۷ بهینه سازی سازه کره‌ای تحت بار گسترده (الف) نمای جانبی (ب) نمای از بالا 
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به روز حسنی- حسین قاسم نژاد مفری- سید مهدی توکلی 
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شکل ۱۸ نمودار همگرایی بهینه‌سازی همزمان سازه کره‌ای تحت بار گسترده 


نتایج 
در این مقاله از روش بهینه‌سازی همزمان برای طراحصی 
شکل صفحه میانی پوسته. به همراه پیکره لاه سخت 
کننده آن جهت رسیدن به سخت‌ترین سازه ساخته شده 
با مقدار مشخصی مصالح لایه سخت کننده, استفاده شده 
است. در این روش در هر گام بهینه‌سازی, بعد از طراحی 
شکل صفحه میانی. بهینه‌سازی لایه سخت کننده انجام 
شده است. برای تولید شکل از فن تولید سطوح نربز 
ستفاده شده است و نقاط کنترلی سطوح به عنوان 
متغیرهای طراحی در مرحله بهینه‌سازی شکل مورد 
ستفاده قرار گرفته‌اند. در بهینه‌سازی پیکره لایه سخت 
کننده از مدل مصالح ساختگی همسانگرد با اعمال تاوان 
ستفاده شده و چگالی مصالح هر یک از اجزا محدود که 
خود تابع پارامترهای هندسی مواد ریزسلولی تشکیل 
دهنده مصالح هستند. به عنوان متغیر طراحی در نظر 
گرفته شده‌اند. از روش مجانب‌های پویا برای حل مسئله 


بهینه‌سازی هر دو زیرمسئله شکل و پیکره استفاده شده 


راخ 


ست. 

درگذشته به خاطر زیاد بودن تعداد متغیرهای 
طراحی در مسائل بهینه‌سازی پیکره از روش‌های معیار 
بهینگی برای حل مسئله بهینه‌سازی استفاده می‌شد. در 
سال‌های اخیر استفاده از روش‌های برنامه‌ریزی ریاضی 
کن انش شیک ان انا بیان مارل له انسخرفر اس 
مقاله نیز از یکی از روش‌های برنامه‌ریزی ریاضی به نام 
قخاش‌ ها وا باه قاشع همان طوی که مات 
شد با وجود زیاد بودن تعداد متغیرهای طراحی این 
روش قادر به حل مسئله می‌باشد. 

در اين مقاله بهینه‌سازی همزمان شکل و پیکره 
برای رابطه سازی حداکثرکردن سختی سازه درنظر گرفته 
شده است. به‌عنوان کارهای آینده می‌توان معیارهای 
طراحی متفاوتی را با توابع هدف و قید دیگر درنظر 
گرفت. به عنوان مثال می‌توان مسائل کمانش, اعمال قید 
تنش و بهینه‌سازی تحت بارهای دینامیکی و پا تحت 


ترکیب بارها را مورد بررسی فرار داد. 
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